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We discuss analytic approximations to the ground state phase diagram of the homogeneous Jaynes- 
Cummings-Hubbard (JCH) Hamiltonian with general short-range hopping. The JCH model describes e.g. radial 
phonon excitations of a linear chain of ions coupled to an external laser field tuned to the red motional sideband 
with Coulomb mediated hopping or an array of high-Q coupled cavities containing a two-level atom and pho- 
tons. Specifically we consider the cases of a linear array of coupled cavities and a linear ion chain. We derive 
approximate analytic expressions for the boundaries between Mott-insulating and superfluid phases and give 
explicit expressions for the critical value of the hopping amplitude within the different approximation schemes. 
In the case of an array of cavities, which is represented by the standard JCH model we compare both approxi- 
mations to numerical data from density-matrix renormalization group (DMRG) calculations. 
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In recent years there is a growing interest in quantum op- 
tics systems as experimental testing ground of fundamental 
models for quantum many body physics and quantum simu- 
lation. The most prominent example are certainly ultra-cold 
atoms in optical lattices which are almost ideal repre- 

sentations of various types of Hubbard models H 0> S la ED ■ 
However due to their finite mass, atomic systems represent, 
with few exceptions, only models with explicit particle num- 
ber conservation. On the other hand different quantum op- 
tical systems employing photons or quasi-particles such as 
phonons have been suggested recently, which are not limited 
by this constraint. For example an array of coupled high-Q 
micro-cavities containing a two-level atom and a photon is 
described by the Jaynes-Cummings-Hubbard model (JCHM) 
JH IH Oil [Tj]] . It is a combination of two well-known sys- 
tems, the Jaynes-Cummings model lfl2i [l3ll describing the 
coupling of a single two-level system to a bosonic mode and 
the hard-core Bose-Hubbard model Hill which describes the 
interaction and tunneling of bosons on a lattice. Recently 
we have shown that a modification of the JCHM can also be 
implemented in a linear ion trap, which has the advantage 
of an easier experimental realization since the required pa- 
rameter regime is already realizable with current technology 
lfl5ll . A large variety of analytic and numeric methods was 
applied to the JCHM and related models, providing profound 
results for the phase diagram and other ground state quantities 
|H d d El E3 EH El HI . In the present paper we show 
that in the strong interaction limit and near commensurate fill- 
ing simple approximate analytic solutions of the JCHM can be 
found if there is translational invariance, i.e. for an infinite ho- 
mogeneous system or periodic boundary conditions realizable 
e.g. with ions in a race-track Paul trap design. These solutions 
provide a good analytic approximation to the full ground state 



phase diagram. 

This paper is structured as followed. In section Q] we 
will briefly summarize the main properties of the Jaynes- 
Cummings model, together with other important quantities 
needed later on. In section [TT] we introduce two different ap- 
proximation schemes, both giving analytic results for the crit- 
ical hopping amplitude for the Mott insulator (MI) to super- 
fluid (SF) transition. In section [III] we apply both approxima- 
tion to the simple cubic nearest-neighbor JCH model describ- 
ing an array of coupled cavities and to the special case of an 
linear ion chain. 



I. THE JCH MODEL 

In this section, we will shortly review the main features of 
the JCH model defined by the Hamiltonian 



+ a j+d a 3 



(1) 
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and discuss the main quantities needed in order to calculate 
its phase diagram. The system ([T]i comprises bosonic and 
spin degrees of freedom, the specific interpretation of which 
depends on the actual physical system. Depending on the 
implementation, and a,j describe the creation and annihi- 
lation of a photon (phonon) at the jth cavity (ion), are 
the spin flip operators between the internal states of the atom 
(ion) and A is the transition energy of the atom (the detuning 
of the external laser field from the red motional sideband). 
g describes the cavity-mediated atom-photon coupling (the 
phonon-ion coupling in the Lamb Dicke limit) and u) is the 
cavity resonance (the local oscillation) frequency. Between 
separated cavities (ions), there is a photon (phonon) transfer 
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described in ([T]) by the distance dependent hopping amplitude 

In the limit of vanishing hopping td = 0, the resulting 
Jaynes-Cummings model can easily be diagonalized. In this 
case, all sites j decouple and become independent. Since the 
total number of excitations Nj = dz-hj+a^ a J on every site j 
is a constant of motion, the local JC Hamiltonian block diago- 
nalizes. Within each two-dimensional excitation subspace, the 
eigenstates can easily be found. Adapting the notation from 
flal . these are given by 



\±,n) 



[x„T(w-A)]|t,n-l)±2 5 ^|l,7i) 



a r t|T,n-l)±/3±||,n) 



(2) 
(3) 



with Xn = \fW 



Ang 2 and n > 0, and the eigenen- 



ergies are 



E 



A 



-U) 1 



(4) 



For n = 0, the ground state is non-degenerate and 
given by | — ,0) = ||,0) with E = 0. Here, the state 
|t,n—l) describes an atomic excitation together with 
n — 1 bosonic excitations; \[,n) is the state with the 
atom in the ground state and n bosonic excitations. In 
the strong interaction limit g 3> |A — u\, the energy gap 
AE n = E£ — E~ = x n ~ 2g^/n is large compared to any 
other energy scale in the system and thus the excited states 
\ + ,n) do not contribute to the ground state. 

For the following discussion it will be useful to consider the 
action of a single bosonic creation or annihilation operator on 
a given JC eigenstate |±, n). Defining 



fa uiPn+i ± y/n+l (3 n z a n+1 n > 
a7 n = 



(5) 



B± = \^ i a nK+i T y/n+1 K <+i n > (6) 



n = 



n < 1 





± _ J Vn - 1 atPn-i T >/n ^ n > 1 
±^^,i n < 1 



(7) 



, (8) 



the action of and a on the state |±, n) can be seen to be 

a) \±, n) = At | + , n + 1) + B± |-, n + 1) (9) 
a \±, n) = Ct |+, n - 1) + D± |- n - 1) , (10) 



i. e. a> and a connect the manifold of states |±,n) to the 
manifolds |±, n + 1) and |±, n — 1) respectively as expected. 

In order to calculate the phase boundaries of the Mott insu- 
lating lobes for the JCH model, we will follow the usual route. 



Since the total number of excitations in the system 



N 



(id 



commutes with the full Hamiltonian ([]]), it is enough to treat 
the system for a fixed number of excitations. The boundary of 
the 7ith Mott lobe can then be determined by calculating the 
total energy E(N) for TV = nL- 1, N = nL and N = nL+ 1 
excitations in a system with L sites. The chemical potential 
then reads 



(it = ± EinL ± 1) - E(nL) 



(12) 



where the plus sign belongs to the upper boundary of the Mott 
lobe and the minus sign to the lower one. For td = , fi^ can 
be calculated straight forwardly. Starting with the energy for 
N = nL excitations with n being an integer, i.e. for a com- 
mensurate number of excitations, it can be seen that due to the 
nonlinear dependence of the single-site energy E~ on n, the 
excitations will distribute equally over the whole lattice. The 
ground state is therefore given by n — {n, n, . . . , n}. Now, 
when adding (removing) a single excitation from the whole 
system, the ground state is given by {n ± 1, n, . . . , n}, where 
we have ignored the degeneracy of the state since we are only 
interested in the energy and the system is homogeneous. With 
this, the energies at td = can be written as 



E(nL-l) = (L-l)E- + E-_ 1 

E(nL) = LE~ 
E{nL + l) = (L - 1)E~ 



and the chemical potentials evaluate to 

M« = E n+l — E n 

= u — + (1 - o n0 )— + d nfi — - — , 



for any n and 



E n ~ E n-l 



(13) 
(14) 
(15) 



(16) 
(17) 

(18) 
(19) 



for 7i > 0. Thus for a commensurate number of excitations 
the system displays particle-hole gaps. Since fj,~ +1 = fJ,£ 
the chemical potential for non-commensurate total number of 
excitations between N = nL and N = (n + i)L is the same, 
corresponding to a critical point. For non-vanishing tunneling 
the critical points extend to critical regions. 

The simplest numerical method to obtain a qualitative phase 
diagram is the so called mean field theory. As described for 
instance in fiol [Til l2ll l23ll . mean field theory can be imple- 
mented by introducing an order parameter ^, which in our 
case is chosen to be homogeneous and real valued. Decou- 
pling the hopping term by using 



*(*] 



6,-1*1 



(20) 
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the whole JCH Hamiltonian (fl~|i in the grand-canonical ensem- 
ble uncouples in real space with a local Hamiltonian being 

H MF = (w-(jt) a) a + (A - n) a + a~ + g (a^" + aa + ) 
-2J*(a f +a) +2J|*| 2 . (21) 

At this point, we omitted the spatial index because the prob- 
lem is purely local. The modified hopping amplitude J = 
— J2d td gi ves the effective coupling within the mean field 
scheme. The phase diagram is now found by diagonalizing 
the mean field Hamiltonian d2~TT i either exactly by means of 
perturbation theory or numerically, setting an upper bound for 
the maximal number of bosonic excitations in the system. The 
ground state energy is then given by min E[^>] and the MI is 

distinguished from the SF by the value of vf" for the minimal 
energy. For ^ = 0, the system is in a MI state, for > 0, 
the ground state is superfluid. This sets the point of the MI to 
SF transition. It should be mentioned at this point, that this 
method gives inadequate results in one dimension (D = 1) 
but is exact for D — ► oo. Additionally, the effective hopping 
J must be larger than zero to yield useful results. 



II. APPROXIMATIVE DETERMINATION OF THE PHASE 
BOUNDARIES 

A. Effective strong-coupling model 

From the discussion above, it can be seen that the phase 
boundaries are defined by the closure of the particle-hole gap. 
In the present subsection, we will derive effective Hamilto- 
nians in the strong-coupling limit for the calculation of the 
upper and lower chemical potential of the nth Mott lobe, 
allowing to calculate the particle-hole gap in first order of 
the hopping amplitudes td- To do so, we employ degenerate 
perturbation theory using Kato's expansion as summarized 
in ll24ll in first order with H c h = WV . This procedure is 
equivalent to the polariton mapping considered in ITU l22ll . 
First, we note that according to eq. (|4]i the state |+,n) is 
separated by a large energy gap from the ground state |— , n). 
Thus, |+, n) can be completely neglected in the following as 
already mentioned in lflll[l6tl . 

When looking for the energy of the ground state with N = 
nL, from perturbation theory, no 1st order contributions are 
present. So, the Hilbert space per site is one dimensional, 
consisting of the single state |— ,n). Thus, up to first order, 
the energy is given by E{nL) = LE~ . When adding an ex- 
citation, the local Hilbert space increases; now, (locally), the 
two states | — , n) and | — , n + 1) need to be taken into account. 
So, in this limit, the system for an additional particle can be 
understood as a system consisting of effective spin | parti- 
cles. We will identify the states |f|~) with the state |— , n + 1) 
and |4) with |— ,n). In order to derive the effective spin | 
model, one has to look on the action of the hopping operator 
aj +1 cjj on the states in the Hilbert space. Using equations (O 
to (O and neglecting the contributions from the states | + , n) 



and |+, n + 1), the hopping operator d^ +1 dj acts as 



a j+ i a 3 



Mi = B n D n+l 



'3+1 



(22) 



within the considered subspace. Therefore, by introducing 
spin operators af the hopping term is equivalent to a near- 



est neighbor spin-spin interaction with 



a j+l a 3 



B n D n+1 aj +1 a j 



(23) 



Together with the energy of the system, one can thus write an 
effective Hamiltonian describing the upper boundary of the 
nth Mott lobe 



H = E nY,°J°t +E n+lY,*Pi 

3 3 



a 3 :(7 3+d 



(24) 



d j 

This Hamiltonian is equivalent to 



H = (L - 1)E~ 



E 



n+l 



U 3 a 3+d 



(25) 



since we are at fixed magnetization with only one spin point- 
ing upwards. This Hamiltonian can be further simplified, by 
using a Jordan- Wigner transformation mapping the spin op- 
erators aj onto fermionic operators Cj and subsequently per- 
forming a Fourier transformation 



Cfe. 



(26) 



Then, the ground state wave function factorizes, since the 
Hamiltonian decouples in momentum space 



H = (L — 1)E„ 



J n+l 



kd 

2B-D- +1 ^ d ^cos(27T— )c[c fc . (27) 



This model is equivalent to free fermionic particles with hop- 
ping amplitudes given by td- In momentum space, a single 
fermion will occupy the mode with lowest energy. Thus the 
total energy of the single particle, and therefore the total en- 
ergy of an additional excitation on top of the nth Mott insula- 
tor in the JCH model is given by 



E(nL + 1) 
where 



(L 



De: 



E 



n+l 



F n {k'), (28) 



F n (k) = 2 B~D- +1 ^ cos(2^— 



(29) 



and the momentum mode k' is chosen such that F n (k)L, is 
minimal. It should be mentioned that the product B~D~ +1 is 
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positive for any (A, oj, n), so the momentum mode is purely 
determined by the minimum of ^2 d t& cos(27r-^). 

To calculate the energy for a hole in the nth Mott insula- 
tor, we follow exactly the same route. Now, the state |JJ.) is 
associated with |— ,n — 1) and |f|~) with |— ,n). The hopping 
operators act as 



and the effective Hamiltonian is given by 



(30) 



h = e-_ x J2 + e- J2 °P7 

i i 

+ B-_ t D- J> £ (°U°7 + *t*7+*) ■ w 



Here the magnetization consist of one spin pointing down- 
ward. Again, after making use of a Jordan-Wigner transfor- 
mation and subsequently a Fourier transformation, the energy 
of a single hole is given by 

E(nL - 1) = (L - 1)E~ + + F n _i(fc"). ( 32 ) 

where the same condition holds for k". Now, putting the cal- 
culated energies d28l ) and d32t together, the chemical poten- 
tials and therefore the boundaries of the nth Mott insulating 
lobe can easily be derived. They are given by 

/4 = E n+i - K + 2B-D- +1 £ td cos(27r^) (33) 

d 

M« = E~ - E~_ x - 2B-_ ± D- U cos(27r^). (34) 

d 

where k' (k") is chosen such that /j+ (fc') (/i~ (k")) is minimal 
(maximal). This result generalizes the findings from lfl6ll22Tl 
to arbitrary hoppings td- 



B. Fermion approximation 

In this subsection, we will apply an even simpler, but 
not that obvious approximation. When looking at the JCH 
Hamiltonian (fl}, it can be seen that all terms are quadratic. 
These kinds of models are in general suited for an exact 
solution by means of a Fourier transform. The problem 
at this point is however, that the commutation relations of 
spin operator cr^ are not as simple as that of bosons or 
fermions. This limits the applicability of a Fourier transform, 
since the operators in momentum space will not obey the 
same commutation relation as in real space |26|. The usual 
step of a prior Jordan-Wigner transformation, transforming 
the spin operators to proper fermionic operators, is not 
applicable in this case, since the interaction part is linear 
in the spin operators, so the Jordan-Wigner factors do not 
cancel out. Thus both transformations cannot be carried out 
exactly without increasing the descriptional complexity of the 



problem. Nevertheless the Hamiltonian can be diagonalized 
by a Fourier-transform in an approximate way. 

As said above, all modes decouple at td = 0. For this rea- 
son, the spin-operators are in this limit equivalent to fermionic 
operators. If we assume that this replacement also holds for 
small values of td, the JCH model (fl} can be rewritten in a 
fermionic approximation 



(35) 



Here the spin operators a + (&~) are replaced by fermionic 
operators & (c). Within this approximation, a Fourier trans- 
form of both, the bosonic and fermionic degrees of freedom 
can be easily accomplished via 



Vl 



Ck- 



(36) 
(37) 



Here and are operators in momentum space. Doing so, 
the JCH Hamiltonian transforms to that of uncoupled Jaynes- 
Cummings systems 



H = > ujka' k a,k 
k 



9 



E(4 



a\c k 



with 



kd. 



LU k = OJ + 2 ^2 td C0S(27T— ) 



(38) 



(39) 



The ground state in any mode is given by the Jaynes- 
Cummings ground state © with frequency u>k. The energy 
of mode k with n excitations is 



Eu = (1-<W 



nujk 



(40) 

Since the total number of excitations in the system 



N 



a\a k 



c\c k 



(41) 



commutes with the Hamiltonian d38l , a common basis can 
be chosen. Thus the full solution of d38l for a fixed total 
number of excitations N = nL is given by the distribution 
?i = {n kl , n k2 , . . . } of N excitations on L momentum modes 
with minimal energy En[ti\ = J2 k E k k together with the 
constraint J2 k n k = N. Note that the number of momentum 
modes L is equal to the number of sites. 

When constructing the phase diagram, the energy of N = 
nL — 1, N = nL and N — nL + 1 excitations needs to 
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be calculated. In the limit of vanishing hopping (t = 0) 
and for commensurate filling, i.e. N = nL, the distribution 
of occupation numbers which has the lowest energy is again 
?i = {n, n, . . . , n}. This corresponds to a MI state with an 
integer number of excitations on every lattice sites. The phase 
is gapped with a particle-hole gap as described in section U 
When t is increased the ground state remains the same, but 
the gap closes and a quantum phase transition occurs from the 
MI to the SF phase at some critical value of t. The only re- 
maining thing in order to calculate the chemical potentials is 
to find the momentum mode where the addition (removal) of 
an excitation gives the maximum (minimum) reduction (in- 
crease) in the total energy. This yields 



(42) 

^-=El-El~\ (43) 

where kl (k) is chosen such that (u~(k)) is minimal 

(maximal). The actual values of k and k' depend mainly on 
the sign of the hopping amplitudes ta- 



rn. APPLICATION TO SPECIFIC REALIZATIONS OF 
THE JCH MODEL 

After having introduced the two approaches used in this pa- 
per, we will apply them to the case of the simple JCHM with 
positive effective mass and nearest neighbor hopping and to a 
modified model describing the physics of a linear ion chain. 
The case of the simple JCHM essentially serves as a testing 
ground for our approximation schemes including a compari- 
son of the analytic results to numerical data from DMRG and 
mean field calculations. Later on, the generalized JCHM will 
be treated by both approximations giving analytic results for 
the phase diagram in a wide range of parameters. 



A. JCHM with positive effective mass and nearest-neighbor 
hopping 

Without loss of generality, we will specialize here on the 
case discussed in 112011 . The Hamiltonian of the JCHM in this 
case is given by 



H = uj ojfij + A J2 + 9 J2 + 

3 3 3 

-*Z)( a 3 a J - + i+a 5+i a J - )- ( 44 > 



with uj = A. Comparing with the Hamiltonian eq. (Q]i, one 
notes that the hopping amplitudes satisfy td = —tSdx- 

For the calculation of the chemical potentials, we first have 
to determine the momentum modes k' and k" which con- 
tribute to the energy. For u> = A, the coefficients in (O are 
&n = ~^2 = Pn ar, d therefore 



B~ = 




n > 
n = 



D 



n+l- 



(45) 



^—effective model 
■ - -fermion approximation 
x DMRG 
mean field 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

hopping t/g 



FIG. 1: (Color Online) Comparison of ground state phase diagram 
of the ID JCHM J® obtained by DMRG (red crosses, from Rossini 
and Fazio PRL 99 186401 (2008)) as well as mean field results (dot- 
dashed line) with the prediction from our approaches (solid line: 
strong-coupling effective Hamiltonian; dashed line: fermion approx- 
imation) for A = uj — 1 and g = 1. Taking into account the 
simplicity of both approaches, the agreement with the DMRG data 
is rather good while the mean field predictions are rather poor as ex- 
pected for ID systems. The critical hopping amplitudes estimated 
from the DMRG data agree surprisingly well with those predicted 
within the fermion approximation, although the shape of the Mott 
lobe is different. 



With this, the function F n (k) is given by 



* — O n .O 



L' 



The both chemical potentials have its minimum (maximum) 
at k = 0. Putting everything together, the phase boundaries 
of the nth Mott lobe, calculated using the effective strong- 
coupling model read 



1 



for any n and 



1 , l-<Jnl 

M„ = <•> - + 2 — 



-Xn ~ t 



2-5. 



(47) 



Xn-l + t 



(x/n + y/n-l)' 



2 - S nl 



(48) 



for n > 0. This allows for the determination of the critical 
hopping amplitude < crit where //+=//", which is given by 



tciit/9 — 2 



2y/n - Vn + 1 - y/n-1 



(y/E + Vn + I) 2 + {Vn + Sni + Vn - l) 2 



(49) 

Secondly we apply the second approximation to this model. 
With the given system parameters, the momentum dependent 
phonon energies d39l are given by 



LUk = uj — 2tcos(27r— ) 



(50) 
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and the energy in the fcth momentum mode for a given filling 
n reads (see eq. (|40l i) 



boundaries within these approximation giving explicit ana- 
lytic results for them. 



= (1 - <W 



k k 
nu> — 2nt cos(27r— ) + t cos(27r— ) 
L L 



t 2 cos 2 (27r— ) + ng 2 



(51) 



Finally, following eqns. (l42l and d43l l. the momentum modes 
k' (k) which minimize (maximize) the chemical potentials 
need to be found. In the present case (ti < 0), these are 



k' = and k 



Thus the resulting chemical potentials are 



i+ - u = -2t + tS„o - \A 2 + («+!).( 



+ (l-S nQ )^t 2 + ng 2 , (52) 



for any n and 



- u> = 2t - t8 nl - \[t 2 + 



+ (l-5 nl )^t 2 + {n-l)g 2 , (53) 

for n > 0. A closed form for the critical hopping can be 
found, but is rather lengthy and will therefore be skipped. 

We now compare our analytic results to various numeri- 
cal calculations. Figure[T]shows both analytic approximations 
along with numerical data from DMRG H2Cj| and mean field 
ll23ll calculations, where the modified hopping amplitude in 
the mean field Hamiltonian (|2~TT i evaluates as J = t. From 
the figure, it can be seen that the effective model gives a 
much better agreement with the numerical DMRG data, es- 
pecially the slopes of the lobes agree perfectly at small hop- 
ping. The fermion approximation overestimates the size of the 
Mott lobe. In particular while the lower boundaries are rather 
well reproduced the upper boundaries have the wrong slope. 
Surprisingly though the critical hopping amplitudes seem to 
agree better with the DMRG data than the results obtained 
from the effective strong coupling Hamiltonians. Although 
the fermion approximation is quantitatively worse than the ef- 
fective strong-coupling Hamiltonians, it provides a simple ap- 
proximative solution to the JCHM beyond the mean field level 
which has the advantage of giving a closed form of the ground 
state. 



B. Linear ion chain 

As a second example we consider a linear string of ions in 
an ion trap lfl5ll . where the ions are coupled to an external 
laser field and interact with each other due to Coulomb 
repulsion via phonon exchange. This system is well described 
by a modified JCH model with a specific short range hopping 
with negative effective mass and site dependent parameters. 
First we will shortly introduce the model and give a derivation 
of the corresponding homogeneous limit. Afterwards we 
will apply the both approximations and discuss the phase 



As shown in I15H . the Hamiltonian of a linear string of L 
ions simultaneously irradiated by a laser which is tuned close 
to the red radial motional sideband and in the Lamb-Dicke 
regime is given by 



L-l 



3=0 



L-3 



+ ^ \ & ]+d h 3 + a]a j+d J ■ (54) 

j d=i 

Here a\ and a, describe the creation and annihilation of a lo- 

3 j 

cal phonon at the jth site (ion), are the spin flip operators 
between the internal states of the ion, A is the detuning of 
the external laser field from the red motional sideband, g de- 
scribes the phonon-ion coupling in the Lamb Dicke limit; for 
a precise definition see lfl5ll . The local oscillation frequencies 
<jjj and the hopping amplitudes tj.j+rf are determined by the 
longitudinal and transversal trap frequencies lo z and oj x via 



2uj x 



L-l 

E 

1=0 



1 



1 



Ul\ 



3,3+d 



2oj x 



\Uj - u j+d \ 



(55) 

where Uj are the equilibrium positions of the ions 12511 . For 
sufficiently large L, the equilibrium positions of the ions at the 
center are approximately equidistant, giving Uj = j u, with u 
being the distance of two adjacent ions. 

Let us now discuss the limit of a homogeneous chain ne- 
glecting any boundary effect. In this limit eqs. d55l ) can be 
rewritten for L — > 00, yielding position independent phonon 
energies ujj = —uj and hopping amplitudes tj j+a = td 



1 



2io x v? d 3 
,2 



1 



u,'Z 



2ui x v? 



C(3) = 2fC(3), 



(56) 
(57) 



where t — 2t ^ acts as a small parameter and u> > 0. ((x) 
is the Riemann ^-function. 

One notices a negative oscillator energy — uj and a neg- 
ative effective mass, which is a result of the positive hop- 
ping strength t. This negative mass is the reason, why the 
application of the mean field theory is not that straight for- 
ward. When simply calculating the modified hopping ampli- 
tude J = —t J2 d = ~~ *C(3). the hopping gets negative and 
therefore mean field theory inapplicable. This problem can be 
overcome by first applying a canonical transformation to all 
used operators. The transformation 



(-i)'a 3 - 



(58) 



for the annihilation operator and accordingly to all the other 
operators v-ppf, maps the JCH model (Q]) again to the JCH 
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* 1st order fermion approximation mean field 




0.001 0.01 0.2 0.001 0.01 0.2 0.001 0.01 0.2 

hopping t/g 



FIG. 2: (Color Online) Phase diagram of the JCH model for a linear ion chain for three depicted values A/g = —0.8, 0, 0.8. Shown are the 
upper boundary of the zero filling lobe (always lowest line) and the boundaries of the lobes with filling from 1 to 5 on a double logarithmic 
scale. Beside the used approximations (solid line: fermion approximation, crosses: 1st order effective theory) the results from the mean field 
theory (dot-dashed line) after the canonical transformation are shown. It can be seen, that the fermionic approximation again overestimates 
the phase boundary (compared to the more reliable effective strong coupling theory) but gives a better agreement compared to the mean field 
theory (mind the logarithmic scale). 



model, but with i<j i— > (— l) d td- After this transformation, the 

modified hopping evaluates to J = —t^2 d y d $ = 3t£(3)/4 
being positive. Now the application of the mean field theory 
is straight forward, following the usual route. 

After having introduced the homogeneous limit of the 
model, the approximations introduced in section [XT] will both 
be applied. Starting with the effective strong-coupling theory 
from section III Al the chemical potentials for the upper and 
lower boundary of the lobes are given by eqns. d33l l and d34i >. 
The proper momentum modes k', k" , which minimize (maxi- 
mize) the chemical potentials are both found to be k = L/2. 
This results from the negative mass. Due to the complexity 
of the problem, especially the analytic form of B~ and D~, 
analytic representation of the chemical potentials are left out 
here. They can be found straight forwardly just as in the case 
of the simple cubic JCH model. 

When following the approximative method from section 
III Bl the Hamiltonian for the uncoupled JC models is given 
by eq. d38l l. with the phonon energies being 



2t 



x ^ cos(2tt^) 



d 



d 3 



(59) 



negative. Using the polylogarithm Li„ (x) = J^dLi W one 
can write them in the explicit form 



t Li, 



,2™-, 



Li 



3 e 



-27ri-, 



2C(3) 



(60) 



The minimum value of uj^ 



-7tC(3)/2 is attained for 



k = j as expected from the positive sign of the hopping 
term. The energies for each momentum mode are given by 
the solution (l40b of the JC model and the corresponding 
spectrum is shown in fig. [3] 

From the knowledge of the dispersion relation for different 
fillings, it is now easy to construct the phase diagram. As dis- 
cussed in section UTBl the flat dispersion for t = leads to the 
ground state having an equal number of excitations in every 
momentum mode k. The chemical potentials for t > are 
then determined by the k' and k values, minimizing or max- 
imizing eqns. d42l) and d43l >. When looking at the dispersion 
in figure|3] one recognizes that this is given for k' = L/2 and 
k = 0. So, the chemical potentials are given by 



= E" r +1 - E r } 



a 



EH 



(61) 
(62) 



according to < f39b - Note that since ui = 2t£(3), all t^'s are 



and when using the analytic form, eqs. (l40l > and d60l i. the 
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momentum k/L 




detuning A/g 



FIG. 3: (Color Online) Energies of the JCH Hamiltonian for fixed 
filling n as function of momentum k. Shown are the energies from 
eq. 10 for the five lowest fillings ... 4 (from top to bottom) for 
A = and g = 1, Solid lines: t/g — 0.02; Dashed lines: t/g = 
0.2. One clearly recognizes the minimum at k = L/2 and the flat 
dispersion for t/g — > 0. 
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""9 4/g 



hopping t/g 



FIG. 4: (Color Online) Phase diagram of the JCH model for a linear 
ion chain from the fermion approximation. Boundaries of the Mott- 
insulating lobes (from bottom to top) for n = 2,3, 4. The lobes 
n = and n = 1 are not displayed since they are unbound for 

A -» -oo. 



FIG. 5: (Color Online) Critical hopping amplitude t„ Tlt (A) giving 
the point where the MI to superfluid transition takes place. From top 
to bottom: n = 1 ... 8, all for g=0.05 



typical lobe structure of the Mi-phases with a closing of the 
lobes at some value ^ rlt (A). Whilst the mean field results 
underestimates the extent of the MI regions, our fermionic 
approach overestimates them but with a better agreement 
with the first-order effective strong coupling model compared 
to the mean field solution. The main advantage of the 
fermionic approximation is the easy closed form for the 
chemical potentials as well as the for the ground state and a 
more reasonable agreement of the critical hopping amplitude 
t$j rlt (A) as can be seen from the figure. Figure H] shows the 
full phase diagram of the model as a function of the detuning 
A obtained from the fermionic approximation only. 

The critical hopping amplitude i™'(A) can easily be 
calculated from the analytic expressions for the chemical 
potential given above. Figure [5] shows the dependence of 
the critical hopping amplitude from the detuning A for the 
different MI lobes. One recognizes the unboundness of the 
first lobe, i.e. t^ nt (A) — » oo as A — * — oo. 



phase boundaries of the nth Mott lobe read 

^ 1 



( l((3)t + A 
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1 + 5 n 



-C(3)t + S n0 A 



W + (^C(3)i + A) 2 , (63) 



-IVW + A2 + ^iA. (64) 

Figure|2]shows the resulting phase diagram for three values 
of A comparing the different approaches. One recognizes the 



IV. SUMMARY 

In summary we have presented two simple analytic ap- 
proximations to the phase diagram of the Jaynes-Cummings- 
Hubbard. The first approximation describes the particle-hole 
excitations in the vicinity of the Mott-insulator to superfluid 
transition for a specific filling by a simple effective spin 
model which generalizes the know results to arbitrary short 
range hopping. The second approximation treats the spins as 
fermions which allows for a simple solution of the model by 
means of a Fourier transformation. A comparison of both 
methods to DMRG and mean field data shows reasonable 
agreement to the numerics. The approximative description 
by effective strong-coupling Hamiltonians makes very good 
quantitative predictions for the phase boundaries of the Mott- 
insulating lobes for small hopping and can be straight for- 
wardly written down up to 2nd order. The fermion approxi- 
mation performs also very well for the lower boundaries but is 
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less accurate for the upper ones. It does make however rather 
good predictions for the critical hopping at commensurate fill- 
ings and has the advantage of giving a closed form for the 
ground state in the whole parameter regime. Altogether both 
methods provide quite reasonable results for the phase bound- 
aries compared to numerical results from DMRG simulations. 
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